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ABSTRACT 

We present a high precision, large sample luminosity function (LF) for the Galactic 
globular cluster M3. With a combination of ground based and Hubble Space Telescope 
data we cover the entire radial extent of the cluster. The observed LF is well fit by 
canonical standard stellar models from the red giant branch (RGB) tip to below the main 
sequence turnoff point. Specifically, neither the RGB LF-bump nor subgiant branch LF 
indicate any breakdown in the standard models. On the main sequence we find evidence 
for a flat initial mass function and for mass segregation due to the dynamical evolution 
of the cluster. 

Subject headings: globular clusters: individual (M3) — stars: red giant— —stars: evo- 
lution 
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1. Introduction 

One of the classical tests of stellar evolutionary cal- 
culations is through the comparison of theoretical and 
observed stellar l uminosity functions (LFs) of globu- 
lar clusters (e.g., Renzini fc Fusi Pecci 1988). The 



LF for the stars below the main sequence turn-off 
(MSTO) can be related to the stellar initial mass func- 
tion and be used as a probe of stellar dynamics within 
a cluster. The LF of stars after the MSTO depends 
primarily on the rate of evolution and provides a di- 
rect and straight- forward test of evolutionary calcula- 
tions. Basically the post-MSTO LF measures the de- 
velopment of the hydrogen burning shell and its pro- 
gression through the star. Features in the observed 
LF can be related to interior structure. For example, 
as first noted by Iben (1968), when the H-burning 



shell advances through the composition discontinuity 



left by the maximum penetration of convective envc- 
lope, the luminosity briefly decreases. This leads to 
the so-called LF bump. More subtle observables are 
also present: the slope of the LF below and above the 
LF bump are different because the H-burning shell is 
in the first case passing through a region of varying 
H abundance and later through a region of constant 
H (|Rood fc Crocker 1985|, |Fusi Pecci et al. 1990]}. A 



breakdown in canonica l stellar evolution theory plcn- 



rini fe Fusi Pecci 1988 ) can affect the LF in the region 
of the subgiant branch (SGB).0 Examples include 
extra energy transport b y Weakly Interacting Mas - 
sive Particles (WIMPs) (|Faulkncr fe Swcnson 1993;) 



or the gravitational settling of helium ( Proffit fe Van 
|:lcnBcrg 1991 , Castcllani et al. 1997 , Stranicro et al 



1997), both of which affect the cluster age calibration. 

To explore short evolutionary phases or to search 
for the small changes in the LF produced by reason- 
able modifications to the canonical models requires 
large samples. Indeed, to explore phenomena like 
those cited above, such large samples are required 
that photometry must extend to the crowded in- 
ner parts of clusters. To be useful in constructing 



12 The canonical assumptions include the neglect of phenomena 
like rotation, non-convective matter transport, etc. in stellar 
model calculations. Some of these are certainly present in real 
stars, but they are very difficult to include in calculations. This 
leads to two approaches to dealing with non-canonical phenom- 
ena: one is to try to model the phenomenon ordinarily in a 
simplistic way requiring the introduction of free parameters. 
The other is to make the most stringent possible comparison 
between observations and canonical models — the approach we 
adopt here. 



a LF, photometric samples must be complete and 
accurate — a requirement which has come into conflict 
with the necessity of working near cluster cores. Ad- 
equate samples to identify the LF-bump were first 
obtained for 47 Tuc by King, DaCosta & Demar- 
que (1985) and for 11 other clusters by Fusi Pecci 
et al. (1990) (see Cassisi & Salaris 1997 for a review). 
At lower luminosities LF studies have produced some 
puzzling results. Stetson (1991) in a LF formed from 
combined data from the low metallicity clusters M68, 
M92, and NGC 6397 found an excess of stars on 
the SGB just above the MSTO. Bolte (1994) found 
a similar result for M30 which has been confirmed 
by Bergbusch (1996), Guhathakurta et al. (1998) and 
Sandquist et al. (1999). This excess was especially 
intriguing because that was the "signature" expected 
if WIMP energy transport was important (Faulknci 
fc Swcnson 1993 ). However Sandquist et al. (1996)| in 



the largest sample LF obtained to date did not find a 
similar result for the somewhat more metal rich clus- 
ter M5. 

Clearly larger and better LFs are desirable, and 
these are now possible with the Hubble Space Tele- 
scope (HST) which has the capability to obtain pho- 
tometry near cluster centers. While very large HST 
samples of fainter stars can be obtained, the small 
field of view limits the number of evolved stars which 
can be observed. However, the combination of ground 
based photometry of the outer cluster and HST data 
can yield a large sample LF which spans the range of 
cluster stellar luminosities. In this paper we present 
such an LF obtained for the cluster M3 (NGC 5272). 
This is the largest, most precise, and most complete 
LF yet obtained for a globular cluster. 

2. HST-WFPC2 Data Sample 

2.1. Observations and Data Analysis 

The observations were obtained April 25, 1995 as 
part of the HST program G05496 (P.I. F. Fusi Pecci), 
the first of our programs devoted to testing the the- 
oretical models using luminosity functions of a set of 
Galactic Globular Clusters with different metallicity. 

A set of exposures was secured through UV (F255W 
and F336W) and optical filters (F555W and F814W). 
Preliminary results on UV data has been already 



present ed in Ferraro et al. (1997a) and Laget et al 



1998) ; here we will report on results obtained using 
only the V and / exposures. Total exposure times 
were 400 sec in V (F555W) and 560 sec in / (F814W); 
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each exposure is listed in Table 1 by Ferraro et al. 
(1997a). 

The initial data analysis (including bias, dark and 
flat-field corrections) was made through the standard 
HST pipeline. A master unsharp mask image for 
each filter was made using MIDAS to register all the 
long-exposure images, yielding a median frame with 
cosmic rays statistically eliminated. The search of the 
individual star components was made on the V me- 
dian frame, following the usual procedure with RO- 
MAFOT (Buonanno et al. 1983). Then we performed 
the PSF fitting on each individual frame separately 
using a Moffat function plus a numerical map of the 
residuals to better account for the contribution of the 
stellar PSF wings. The average instrumental magni- 
tudes for each star were transformed to the standard 



Jo hnso n V and / system using the Holtzmann et al 
(19 |95)| recipes. The magnitudes thus obtained were 
eventually shifted to match the independent zero- 
point calibration adopted in Ferraro et al. (1997b, 
hereafter M3CCD). The V,V-I CMD for the HST 
global sample is shown in Figure [ij 

Special care has been devoted to ensure that the 
effects of blends (and crowding) do not significantly 
alter the derived LF. The effects of blending can be 
quite problematic in CMDs of dense regions of glob- 
ular clusters. On the main sequence, blended stars 
can potentially be mistaken for binary stars, while on 
the subgiant branch they be shifted erroneously to 
brighter bins. Blends can mimic incompleteness by 
shifting the two stars off the main loci such that they 
are dropped from the LF. To minimize these effects 
suspected blended objects with multiple stellar com- 
ponents lying off the main branches of the CMD were 
individually examined by eye, exploiting the fully in- 
teractive capability of the package ROMAFOT. In 
this way we were able to check and, if necessary, elim- 
inate spurious objects. The degree to which blends 
contaminate the final sample can be estimated from 
the artificial stars tests described below. Some of the 
artificial stars which were not recovered were lost as 
blends. Thus the degree of incompleteness at a given 
level serves as an upper limit to the blend-fraction. 

2.2. Completeness 

A crucial component in the calculation of an ac- 
curate LF is the determination of corrections for in- 
completeness. In our data crowding is the primary 
source of incompleteness. Moreover, because of the 
large gradient in the spatial density of stars present 



in our frame, the completeness at any level of magni- 
tude must be determined as a function of the position 
with respect to the cluster center. 

In order to quantify the efficiency in detecting stars 
as function of the magnitude and the degree of crowd- 
ing we carried out extensive artificial star tests in the 
V filter frames, following the procedure described be- 
low. 

1 . Each chip was split into two annular regions cen- 
tered on the cluster center: 

INNER: r < 50" including the entire PCI field 
and the innermost part of the three WFC 
fields; 

OUTER: r > 50" covering the outer part of 
the three WFCs. 

M3 is a cluster of moderate central density. 
With a core radius r c = 30" (Djorgovski 1993) 
one would expect the stellar density to be ap- 
proximately constant across the PC as is appar- 
ent from Figure |[ However, the stellar density 
has dropped appreciably at the outer parts of 
the WFs. We have adopted 50" as a bound- 
ary between regions of higher and lower pro- 
jected stellar density. A map of the WFPC2 
field of view is shown Figure |^ with the circle 
(at r ~ 50") delineating the two regions. 

2. A set of artificial stars was generated using the 
Tiny Tim package (Krist 1994), and these were 
added to the original frames in random posi- 
tions. In order to avoid a spurious crowding en- 
hancement, the size of the added sample (about 
300 total artificial images or ~ 80 for each chip) 
was < 10% of the total number of stars origi- 
nally found in each chip. Then the star search- 
ing and fitting procedures were repeated on the 
artifically enriched frames. 

3. An artificial star was considered to be "recov- 
ered" on the basis of the following criteria: 

• the differences in position, Ax, Ay, be- 
tween the input artificial star and the re- 
covered one are Ax, Ay < 1 pixel; 

• the difference in magnitude, Amag, is Amag\ < 
0.3. 

The fraction of artificial stars missed during any 
step of the procedure is a good estimate of the 
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actual number of stars missed at each magni- 
tude level during the reduction procedure. 

4. The completeness factor (<£>) was finally com- 
puted as the numbers of stars recovered (N rec ) 
with respect to the number of stars simulated 
(N a i m ) during each trial. The final value of <& 
has been obtained as an average of 3-6 inde- 
pendent trials for each magnitude bin (0.2 mag 
width) in each of the two regions (INNER and 
OUTER). 

The result of these tests is summarized in Fig- 
ure |a,b. The completeness at the MSTO level (V ~ 
19) is > 90% in both the INNER and OUTER re- 
gions, thus the star counts up to that level have been 
corrected by less than 10%. The corrections for com- 
pleteness were applied separately in each part (inner 
or outer) of each field (PC, WFCs) so that the total 
corrected sample properly accounts for the differing 
degrees of crowding. 

2.3. The Mean Ridge Line of the RGB, SGB, 
and Upper MS 

If there were no photometric errors all of the RGB, 
SGB, and Upper MS stars would lie along a line in 
the CMD. We approximate this line with a mean ridge 
line (MRL) which has been determined by an iterative 
procedure: 

1. A first rough selection of stars belonging to 
the RGB, SGB, MS was performed by eye, ex- 
cluding the horizontal branch (HB), asymptotic 
giant branch (AGB), and blue straggler stars 
(BSS); 

2. This preliminary sample was divided into 0.5 
mag width bins; 

3. For each bin the median V — I color was com- 
puted; 

4. All stars lying more than 5<7 V --/ from the median 
value was rejected. 

Points (3) and (4) of the procedure were iterated 
until the solution was stable. The resulting MRL 
points are listed in Table 1. Only about 2% of the 
sample was rejected in determining the MRL. 

Luminosity functions were constructed for each 
chip by dividing each individual CMD in 0.2 mag bins, 



and simply counting all the stars in each bin within 
±5ov_j of the MRL. 

Figure ^ shows a zoom of luminosity function in 
the Turn-Off region for the total HST sample, before 
and after correction for completeness. Errors have 
been computed accordingly to the following relation 
(see also Bolte 1989): 

.N 1 ' 2 iV • o$. 

-bin = [~ r + (1) 

where N is the number of stars observed in each bin, 
$ is the completeness factor, and o$ is the associ- 
ated error which was determined from the rms of the 
repeated completeness trials. 

3. Matching the Ground-Based Observations: 
the Global LF 

HST data presented here complement the large 
data set we have obtained for M3 during the last 
decades. The photographic sample was originally 
presented by Buonanno et al. (1994) and the CCD 
ground based data have been discussed by Ferraro 
et al. (1993, 1997b). The entire photographic sample 
extends over an area between 2' and 7' from the clus- 
ter center. Here we consider a subsample (hereafter 
referred to as PHOTO) of stars lying in the annulus 
between 3 '5 < r < 6' from the cluster center (region 
Al in Buonanno et al. 1994). In this annulus all stars 
down V ~ 22 have been measured. The CCD data 
cover a field of view 7' x 7' roughly centered on the 
cluster center; stars in the innermost central regions 
(r < 20") have been not measured. In this sample the 
deeper exposures (V and I) only reached to V ~ 20. 
In order to ensure an high level (formally 100%) of 
completeness we include only stars with V < 17 lying 
in the annulus 100"~210" from the cluster center. We 
will refer to this subsample as CCD. The area cov- 
ered by the CCD subsample and the the others fields 
in shown Figure 1 of Ferraro et al. (1997b). 

In the PHOTO and CCD subsets (in analogy to the 
HST sample) we applied the procedure described in 
Section 2.2, selecting bona fide RGB-SGB-MS mem- 
bers taking all stars within ±5er in color from the 
MRL in the V, B-V and V, V-I CMDs respectively. 
The CMDs of all stars selected in the three samples 
(PHOTO, CCD, HST respectively) are shown in Fig- 
ure |5|. For each of these samples we independently 
derive the LF. Star counts have been corrected using 
the respective completeness curve: the PHOTO sam- 
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Table 1 

Mean ridge line for the HST sample. 
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0.079 


15.09 


1.017 


18.86 
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-0.088 


16.94 
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pie has been corrected using Table 7 by Buonanno et 
al. (1994), the CCD sample (for V < 17) is virtually 
complete, thus no correction has been applied, and 
the HST sample has been corrected using the com- 
pleteness curves plotted in Figure [?]a,b. In construct- 
ing the global LF we considered only the magnitude 
range where the completeness was greater than 50%. 
Thus the global LF is restricted to V < 21 by the 
most incomplete region in the HST sample — the inner 
region of the WF4 and WF2 cameras (Figure ||a,b) . 

Since the CCD sample spanned a different magni- 
tude range than the others, we determined a scaling 
factor to produce a proper match. We followed the 
procedure explained in Bolte (1994): let be the 
number of stars, after the proper completeness cor- 
rection, in bin i of the sub-sample k, where k = 1,2,3 
indicates the PHOTO, the global HST, and the CCD 
sub-samples, respectively. Thus, for k = 1,2 the re- 
sulting LF is simply given by the sum 



N 



E 

k=l 



nth 



in the range 12 < V < 21. However, in order to 
include the CCD sub-sample (k — 3) we derived a 
scaling factor s from the range in magnitude (12 < 
V < 17.0) which is common to all sub-samples: 

_ Sfc=l n i,k Si n i,k 
Sfe=l n i,k Si n i,k 

so that the number of stars for each bin in the com- 
bined LF is 

N = s x (AT- +n i>fe=3 ). 

Table 2 lists the final LF. Column 1 lists the mag- 
nitude bins (which are 0.2 mag. in width); column 2 
gives Ni, the total scaled number of objects per 0.2 
mag; column 3 contains the errors in the Ni. Even 
though it extends only 2 mag below the turnoff, this 
LF, including more than 50000 stars, is the most pop- 
ulated LF ever published for a galactic globular clus- 
ter. 

It is useful to compare our LF to that for M5 ob- 
tained by Sandquist et al. (1996, hereafter SBSH), the 
best previous globular cluster LF. Doing this, it is ap- 
propriate to compare sample size in specific parts of 
the CMD. Our "bright" sample which corresponds to 
the upper 4.5 mag of the RGB (or V < 17) contains 



944 stars. For such bright stars SBSH were able to 
work near the cluster center, and their sample in the 
comparable RGB interval was about 1000 stars. For 
fainter stars ground based photometry becomes in- 
creasingly affected by problems of photometric blends 
(see SBSH). Blends can be especially insidious when 
they "move" a star from one part of the LF to an- 
other, because this is just the sort of thing one would 
expect in a breakdown of canonical stellar modeling. 
To avoid blends only the sparsely populated outer 
cluster can be used, and it is especially difficult to 
get adequate samples of the SGB region. Here our 
blcnd-frcc HST sample plays a crucial role. We have 
6085 stars brighter than MSTO; SBSH have 3300. In 
a 0.1 mag bin at turnoff we have 820 stars compared 
to SBSH's 280. 

4. Results 

4.1. Preliminaries 

There are certain issues which must be addressed 
in comparing LFs either with other observational LFs 
or with theoretical LFs. The first is how the samples 
are normalized. Beyond this there is some interplay 
between assumed distance modulus, age, and chemi- 
cal composition. 

4- 1.1. Normalization 

One must determine a "normalization" factor to 
adjust for sample size. In recent LF studies (e.g. 



Stetson 1991] , [Bolte 1994 Sandquist et al. 1996) both 
the lower RGB and the main sequence region have 
been used for normalization. We adopt a different 
approach normalizing to the total number of stars 
brighter than the MSTO at V = 19.10. Operationally 
this has the advantage that the turnoff is well defined 
and because of the large number of stars involved the 
normalization is not much affected by Poisson statis- 
tics. The small imprecision in determining the loca- 
tion of the turnoff (<; 0.1 mag) produces only a very 
small error (<^ 1%) in the normalization. 

4- 1-2. Chemical composition 

For M3 Kraft et al. (1992), found [Fe/H] = -1.47, 
based on detailed fine abundance analysis of high res- 
olution spectra of red giant stars. However, their 
study is not fully self-consistent since while stel- 
lar analysis was performed using model atmospheres 
from the grid of Bell et al. (1976, hereafter BEGN), 
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Table 2 

Luminosity function of the global sample (PHOTO+CCD+HST). 
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2.01 
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2.04 
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0.83 


20.90 


3.59 


2.19 


16.30 
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the reference solar iron abundance was extracted from 
the empirical Holweger and Muller (1974, hereafter 
HM) solar model. This inconsistency was resolved by 
Carretta & Gratton (1997) who use the Kurucz (1993) 
grid of model atmospheres for both solar and stellar 
analysis, giving [Fe/H] = —1.34. The offset with the 
study of Kraft et al. (from which the observational 
material was taken) is simply due to the fact that the 
HM model is ~ 150 K warmer than the BEGN models 
in the line formation region, explaining the 0.13 dex 
difference in the derived [Fe/H]. 

Oxygen in red giants can by affected by interior nu- 
clear processing. Indeed, M3 is one of the several clus- 
ters with an observed anti-correlation between O and 
Na along the RGB (see for instance Kraft et al. 1992). 
This anti-correlation arises because of non-standard 
mixing on the upper RGB. Therefore, the best esti- 
mate for the [O /Fe] ratio is that obtained considering 
only unmixed stars, i.e. giants with [Na/Fe] « 0.0 
dex. Using Kraft et al. (1992) we estimate that 
[O/Fe] = 0.19 ± 0.04 dex based on 5 stars. 

One typically assumes that the other a-elemcnts 
(e.g. Kraft et al. 1993) have the same enhancement 
as O and uses two metallicity parameters, [Fc/H] and 
[a/Fe], in computing models. For many problems one 
can get by with just one metallicity parameter using 
the short cut of Salaris et al. (1993). This makes use of 
the global metallicity ([M/H]) based on the mass frac- 
tion of all elements heavier than helium as compared 
to solar. They found that scaled solar models with 
[Fe/H] = [M/H] mimicked models based on observed 
[Fe/H] with a- enhancement. For M3 [M/H] = -1.2. 
So while ideally we would compare M3 data to models 
with [Fe/H] = —1.34 and [O/Fe] = 0.2, a comparison 
to scaled solar models with [Fe/H] = [M/H] = -1.2 
should be an adequate surrogate. Given typical er- 
rors in metallicity determinations it is appropriate to 
explore a region ±0.3 dex from these values. 

4-1.3. Distance modulus and age 

In principle, fitting cluster main sequences to sub- 
dwarfs with precise Hipparcos parallaxes should lead 
to relatively high precision distance moduli. However, 
M3 was not included in the first round of clusters with 



distances de termined in this way ( Reid 1997 , Gratton 
3t al. 1997). The main reason was that the fiducial 



cluster main sequence ridge line (from Ferraro et al. 
1997b; CCD97) in the V, B - V plane did not reach 
deep enough with the photometric accuracy required 
to apply sub-dwarf fitting to determine distance mod- 



uli. However, using an expanded sample of local sub- 
dwarfs from the newly available Hipparcos catalogue 
(Carretta et al. 1999a) it is possible to give an esti- 
mate for the distance modulus of M3 in the framework 
of the longer distance scale defined by Gratton et al. 
(1997) and Reid (1997). With E(B-V) = 0.02±0.01 
and [Fe/H] = —1.30 (to be consistent with reddening 
and metallicity scales used in Gratton et al. 1997), the 
true distance modulus for M3 is (m — M)q = 15.13 
(Gratton 1998, private communication, adopting the 
Ferraro et al 1997b photometry). A cautionary flag 
to the accuracy of sub-dwarf fitting is set by the fact 
that different groups using this method do not get 
the same result (Pont et al. 1997; Reid 1997; Gratton 
et al. 1997). Small details are important (Carretta 
et al. 1999a), and systematic effects could leave resid- 
ual errors of 0.2-0.3 mag in (m — Al). 

There are other independent estimates of the dis- 
tance to M3. Cudworth (1979) gives a kinematic 
distance for M3 of (m - M) = 14.91. From the 
RR Lyrae Sandage and Cacciari (1990) give values 
of (m — M)o ranging from 14.81 to 15.00 depending 
on which metallicity /absolute magnitude relation is 
adopted for the RR Lyrae. Using the latest metal- 
licity/ absolute magnitude determined from Baade- 
Wesselink distances to field RR Lyrae (Fernley et al. 
1998) and our adopted metallicity ([M/H] = -1.2), 
we get (m-M)o = 14.86 (assuming E(B-V) = 0.01). 
In a recent analysis of the RR Lyrae distance scale, 
Dc Santis (1996) gives (m - M) = 15.03. Marconi 
et al. (1998) find that (m - M) = 14.94 gives a good 
fit to the LF of the lower main sequence. Ferraro 
et al. (1999) have recently obtained (m— M)o = 15.05 
within the framework of a homogeneous re-analysis 
of the evolved sequences of the CMD in a sample of 
about 60 GGCs. They derived the distance modu- 
lus from the comparison of the observed level of the 
ZAHB (Vzahb) an d the theoretical models computed 
by SCL97. Their use of synthetic HBs to determine 
the "observed" Vzahb is a significant improvement 
over earlier applications of the same technique. 

Given the present state of uncertainty we consider 
it reasonable to consider (m — M)v = 14.8-15.2. The 
question of age is intertwined with that of distance. 
The longer distance scales imply smaller ages. For 
example, with the "sub-dwarf distance" given above 
and using stellar models from Straniero et al. 1997 
(private communication), the corresponding absolute 
age for M3 is about 12 Gyr. An decrease in (m — M) 
of 0.07 mag leads to an increase in age of 1 Gyr. Hence 
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ages of 11-15 Gyr for M3 are not beyond question. 

4.2. The Red Giant Branch 

The differential LF for our global sample is shown 
in Figure |[ The dominant determinant of RGB-LF 
is simply energy conservation via the so-call energy 
consumption theorem (Rcnzini & Buzzoni 1986). Ba- 



sically the number of stars in some interval along the 
RGB, Ni, depends inversely on the luminosity at that 
point, Ni oc 1/L. So logiVi oc Mbol, and neglect- 
ing bolometric corrections log Ni should drop linearly 
with magnitude moving up the RGB just as seen. 

Theoretical evolutionary models predict the exis- 
tence of a "feature" in the RGB called the LF-bump. 
This LF-bump marks the evolutionary stage at which 
the H-burning shell passes through the composition 
discontinuity left by the maximum penetration of 



the convective envelope (Iben 1968, Rcnzini & Fusi 



Pecci 1988). As emphasized first by Rood & Crocker 
(1985), the practical detection of such a feature re- 
quires very large RGB samples. Because of this, our 
global sample gives us the possibility to obtain one of 
the most precise determinations of the LF-bump. 

Models also predict that in the integrated LF the 
slope changes on either side of LF-bump. The slope 
change results from the fact that the H-burning shell 
is moving through a region of increasing hydrogen 
abundance below the bump and constant hydrogen 
abundance above the bump. This change in slope 
can more reliably locate the bump than observations 
of the bump itself (Rood & Crocker 1985; Fusi Pecci 
et al. 1990). 

In Figure ^| the LF-bump is very well defined in 
both the differential and integrated LF, and is located 
at V = 15.45 ± 0.05 (in good agreement with the pre- 
vious determinations (Buonanno et al. 1994, Ferraro 
et al. 1997b). Also note in Figure ^|the nomenclature 



we iuse Lu refer Lu the diffeieiiL branches: Olu iiuiueii- 
cla'uii! is Ilia I uf llie llieuiisl rn which llie branches 

are defined in terms of the supposed interior struc- 
ture. There is often confusion, in particular, over the 
term sub-giant branch. To the theorist, and in this 
paper, SGB refers to the transition region between the 
MSTO and the lower RGB where the ridge line in the 
CMD is more horizontal V = 18.2-18.6. Structurally 
this is the branch of thick H-shell burning where the 
degenerate core develops. For historical reasons ob- 
servers often refer to stars along the lower-RGB as 
sub-giants. 



We compare our results to two new sets of theo- 
retical evolutionary tracks/isochrones/LFs which use 
the new OPAL opacities and equation of state (e.g., 
Rogers & Iglesias 1992). This should offer signif- 
icant improvement over earlier models particularly 
for opacity dependent LF features like the LF-bump. 
The first set (the no-diffusion models described in 
Straniero et al. 1997, hereafter SCL97) is available 
for scaled-solar abundances and may be used for any 
given value of a enhancement with the Salaris et al. 
(1993) prescription. The second (Bergbush & Van- 
denBerg 1999, herafter V97) has [a/Fe] = 0.3. The 
V97 models also include improved color-temperature 
relations. On the RGB the V97 models utilize a 
non-Langrangian numerical technique which is not 
optimal for investigations of the RGB bump. This 
technique might be expected to smear and perhaps 
slightly lower the LF-bump. 

In Figure the global differential LF is compared 
to SCL97 and V97 models for (m - M)v = 15.03, 
age = 12 Gyr. For SCL97 models we took [Fe/H] = 
[M/H] = -1.2. The V97 models have [Fe/H] = -1.41 
and [a/Fe] = 0.3 which is equivalent to [M/H] = 
— 1.2. The integrated LF compared to SCL97 is shown 
in Figure ||. The fits are really quite good with the 
region of poorest fit at V ~ 18 just above the SGB. 
That the theoretical and observed LFs agree so well 
for stars brighter than V ~ 18 shows that the hydro- 
gen profile in the region M(r) ~ 0.2-0.4 Mq is given 
fairly accurately by canonical models. 

The values of (m — M), age, and metallicity which 
lead to the best LF fit have the virtues of also agree- 
ing well with the longer Hipparcos distance scale with 
their correspondingly shorter ages and with the abun- 
dances obtained by Carretta & Gratton (1998). Fur- 
thermore, these parameters fit in well with those ob- 
tained from a new re-analysis of the location of the 



LF-bump in a sample of more than 40 GGCs by Fer- 



raro et al. (1999). They find that the previous discrep- 
ancy in the bump location between the observations 
and models is completely removed using new models 
and considering the newer global metallicity scale. 

The good fit to the LF-bump comes at the expense 
of some mismatch at the SGB/lower-RGB transition. 
To explore how varying the parameters affect the fit 
we must quantify the quality of the fits. Formally we 
can perform a Chi-squared test computing 



JV 

E 

i=l 



log(iVobserved)i ~ log(iV mode i) j 



log(A^ 



bscrvcd )i 
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For the best fit to the entire data set (13 < V < 20.5), 
shown in Figure ^, there is a slightly better agreement 
with the SCL97 models (Xg CL9 7 = 0.11 and Xv97 = 
0.13). 

One could conceivably prefer parameters which do 
not give the best overall fit. The factors which de- 
termine the luminosity of the "theoretical" LF-bump 
have been discussed in Fusi Pecci et al. (1990), Chi- 
effi and Gratton (1986) and more recently and exten- 
sively in Cassisi & Salaris (1997). The most impor- 



tant factors are the opacity in the 1-2 x 10 6 K range 
and the degree that convective elements undershoot 
the neutral buoyancy point of the convective enve- 
lope. Elsewhere the post-TO LF depends mostly on 



the rate of fuel consumption (but see §4.3 below), and 
one might expect it to be a considerably more robust 
prediction of the models than the location of the LF- 
bump. Given this one might tolerate some residual 
error in the LF-bump location to achieve a better fit 
elsewhere and, thus, exclude the bump region when 
measuring the quality of the fit. If we do so for the 
fits shown in Figure @ we find xlci^jcbunip = °- 01 
and Xv97_xbum P — 0-02- 

In Figure ^ we explore the effect of varying the dis- 
tance modulus, age and chemical composition. Since 
the quality of the fit is essentially the same for SCL97 
and V97, we show results only for the SCL97 mod- 
els. We have found it impossible to fix parameters so 
that the model fits the SGB region and the LF-bump 
equally well. The best fit in the SGB region results 
from a higher age of 14Gyr. However the good fit 
in the SGB region brings with it a worse fit for the 
LF-bump. The fit to the LF-bump can be restored 
at an age of 14Gyr by adopting a higher metallic- 
ity and (m — M) smaller than we consider plausible 
(Figure |b). The upper end of the allowable (m — M) 
range gives good fit to the LF-bump (Figure ||c) at the 
expense of a slightly worse fit to the SGB region and 
a small age (lOGyr). The Chi-square values both in- 
cluding and excluding the LF-bump region are given 
in Table 3. 

From the ^-values listed in Table 3 we conclude 
that considering both the two main observables of 
the LF (the LF-bump and the SGB-rise) we can find 
a good agreement (with both SCL97 and V97 mod- 
els) when we take age = 12, (m — M)y = 15.03 and 

[M/H] 1.2 (see Figure g). On the other hand, 

if we consider that some residual error still might be 
present in the theoretical prediction of the LF-bump 
location, we get the best agreement between theory 



and observations if age = 14, (m — M)v = 14.80 and 
[M/H] ~ -1.2 (see Figure |a). 

A full discussion of the age should also bring in 
the isochrone fits. These are complicated both by the 
color calibration and the color-temperature relations. 
Probably because of this our isochrone fits using the 
ages and distances determined from the LF give no 
further indication of whether the age of M3 is closer 
to 12 or 14 Gyr. 

4.3. The Lower Red Giant Branch and Sub- 
Giant Branch 

The lower-RGB & SGB LF can be especially in- 
teresting. This is the region in which there has been 
some indication of a breakdown in the canonical mod- 



els for very low metallicity clusters (Stetson 1991 



Bolte 1994 


) and recently confirmed for M30 (|Berg- 


bush 1996, 


Guhathakurta et al. 1998, 


Sandquist et al. 


19991). It 


is the region affected by non-canonical 



assumptions like WIMP energy transport (Faulknci 



k Swcnson 1993, 


DeglTnnocenti, Weiss, & Leone 


(1997) 


), and helium settling 


Proffit & VandenBerg 


1991 




Straniero et al. (1997)" 


, both of which could 



reduce cluster age estimates. Rapid stellar rotation 
(VandenBerg, Larsen, & de Propris 1998) and [a/Fe] 
could also affect the SGB region of the LF. We illus- 
trate all of these in Figure |ll| (see also Figs. 2 & 3 of 
Vandenberg et al. 1998). The predicted changes are 
all small, and age differences can mimic these changes. 

As shown in Figs. ^ & |[ the the theoretical and 
observed LFs agree well in the SGB region and lower- 
RGB. Indeed the fits are pretty good under a wide 
range of assumptions (Figure ||). We should note that 
for the excess of lower- RGB stars in M30 cited above, 
the observed LF lies above the lower-RGB at essen- 
tial every bin. For M3 our "bad" fits correspond to 
only one or two mismatched bins. The fits can be es- 
sentially perfect (Figure ^|a) if we accept some error 
in the bump location. Basically, there is no indica- 
tion for a breakdown in the canonical models. SBSH 
found a similar result for M5. 

On the other hand the observational consequences 
of plausible non-canonical effects are small so even 
our good fit may not place very tight constraints on 
their magnitude. We do not attempt to place such 
constraints now because we feel that a traditional LF 
analysis may not be the appropriate way to analyze 
the SGB. By projecting the star count information 
onto the "magnitude" axis one loses the information 
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Table 3 

Chi-square results between observations and SCL models 



13 < V < 17 



17 < V < 20.5 



13 < V < 20.5 



[M/H] = -1.2 
(m-M) v = 15.03 
age = 12 Gyr 

[M/H] = -1.2 

(m - M) v = 14.80 

age = 14 Gyr 

[M/H] = -0.82 
(to - M) v = 14.65 
age = 14 Gyr 

[M/H] = -1.3 

(to - M) v = 15.20 

age = 10 Gyr 



0.101 



0.194 



0.103 



0.112 



0.014 



0.008 



0.026 



0.018 



0.114 



0.201 



0.130 



0.131 



contained in the color distribution of the stars. In a 
future paper we will report on our efforts to study 
the distribution "along" the SGB (see for example, 



Bcrgbush & VandcnBcrg 1997), and provide quanti- 
tative estimates as to the degree that non-canonical 
assumptions can be ruled out. 

4.4. Below the Main Sequence Turnoff 

Below the MSTO the LF primarily probes the stel- 
lar initial mass function (IMF) and cluster dynamics. 
Having been acquired primarily to study post-MSTO 
evolution, the data we describe in this paper reaches 
only ~ 2 mag below the turnoff which corresponds to 
a mass significantly above the lower mass cutoff for 
main sequence stars. Deeper LFs (based on almost 
25,000 stars) exploring this region are presented in 
Marconi et al. (1998) and Carretta et al. (1999c). 

It is widely recognized that stars in dense systems 
(like GGCs) are subject to energy exchange through 
stellar encounters and rapidly evolve toward a state 
of equipartition of energy (so that stars of lower mass 
will have higher velocities). As a direct consequence 
low mass stars should have larger average distance 
from the center and their distribution will be less cen- 



trally concentrated than that of higher mass stars. 
This mass segregation effect is expected to produce 
an observable effect on the LF if a sufficiently large 
radial region of the cluster is sampled. 

Such mass segregation has already been found in 
GGCs surveyed by HST: NGC 6397 (King et al. 
1995), 47 Tuc (Paresce et al. 1995 ), M15 (De Marchi 
& Paresce 1996), NGC 6752 (Ferraro et al. 1997c). Pi- 
otto et al. (1997) give an extended discussion of this 
topic. Mass segregation has been found in M30 using 
ground based data (Bolte 1989). 

Even with our limited MS coverage, we can see ev- 
idence for mass segregation in M3. In Figure [l^ we 
show a comparison of the LFs in three separate radial 
regions: r < 20"; 50" < r < 100"; 210" < r < 360". 
As above the LFs have been normalized to the total 
number of stars brighter than the MSTO. The inner 
regions show a significant drop of star counts over 
the magnitude range V ~ 20 — 21.5. We interpret 
the observed depletion as due to the mass segrega- 
tion effect. The depletion is obvious even in the small 
interval of masses sampled. The effect shown in Fig- 
ure |l2] cannot due to residual incompleteness of the 
samples. In order to completely remove the effect we 
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would have to assume a level of completeness which 
is far smaller than the estimated value. For exam- 
ple in the innermost sample r < 20" (which is almost 
completely contained in the PC) the measured com- 
pleteness must be almost a factor 2 in error (dropping 
from 0.8 down to 0.45) in order to remove the effect. 

Is the observed decrease in the star counts con- 
sistent with the theoretical expectation? In order to 
derive the relative distribution of stars of different 
mass at different distances from the cluster center we 
have used the formula given by King et al. (1995) fol- 
lowing the procedure as explained in DeMarchi and 
Paresce (1996). The faint end of the LF should be 
depleted in the core region with respect to the outer 
regions by the factor / = (po/Pcxt)™ 2 ^™ 1 " 1 where p 
and p ox t are the stellar densities in the core and in the 
external region respectively, and mi and mi are the 
stellar masses at the faint extreme and TO level of 
the LF, respectively. We use the tables of SCL97 and 
assume (m — M)v — 15.0 to relate observed magni- 
tudes to mass — the MSTO at V = 19.1 corresponds 
to My ~ 4 and a mass of ~ 0.84 M Q , V = 20.5 cor- 
responds to My = 5.5 and a mass of ~ 0.73 M©, 
and V = 21.9 corresponds My = 6.95 and a mass of 
~ 0.61 M©, respectively. 

The spatial densities can be approximated by the 
Eqn. 5 in King et al. (1995). In the case of M3 
the external region lies at r ext ~ 290", and taking 
r c ~ 30" (Djorgovski 1993), we estimate that the ra- 
tio (po/pext) is ~ 900. Adopting the mass values ob- 
tained above, the density of low-mass stars should be 
a factor 900 3 ~ 2.4 less than that in the central. 
This value is compatible with the observed depletion 
which is 975/548 - 1.8. 

We can also explore the stellar initial mass function 
(IMF) using our data. The IMF is normally approx- 
imated by a power-law mass spectrum of the form 
<j)(M)dM = M'^+^dM where <f>(M)dM is the num- 
ber of stars in the mass range M to M + dM. The 
observed slope of the IMF measured in terms of x 
should vary with radius in a way consistent with a 
underlying global slope. To verify this we use the dy- 
namical models of Pryor, Smith, & McClure (1986). 
We compute x using Eq. 3 of Bolte (1989). The two 
magnitude intervals used are 18.55 < V < 20.55 
and 20.55 < V < 21.95. Figure || shows the ob- 
served x plotted as a function of radius. We see 
that a; (observed) varies r in a way consistent with the 
x(global) ~ 0. This is close to the value we inferred 
from our total sample. 



5. Summary 

We have obtained a y-band luminousity function 
for the globular cluster M3 using a combination of 
ground based photographic and CCD observations 
and HST observations of the cluster center. The sam- 
ple is the largest, most complete LF every obtained 
for a globular cluster. It is free from the problems of 
photometric blends which plague ground-based LFs. 
In the crucial MSTO/SGB region our sample is more 
than three times larger than previous LFs. 

Following a traditional LF function analysis we find 
no surprises — canonical stellar models seem to fit the 
data well. In particular, 

• The location of the RGB LF-bump agrees well 
new theoretical models incorporating the OPAL 
opacities and enhancements of the a-elements. 
There is no indication of significant "under- 
shooting" by the convective envelope. The best 
overall fit is achieved with age = 12, (m — 
M) v = 15.03 and global metallicity [M/H] ~ 
-1.2 (equivalent to [Fe/H] = -1.34 and [O/Fe] = 
0.2). On the other hand, if we consider that the 
physics determining the bump location might 
have more uncertainty than that determing SGB 
evolution, we might tolerate a worse fit for the 
bump to achieve a better fit for the SGB. In 
that case we find age = 14, (m — M)v — 14.80 
and [M/H] 1.2 

• The slope of the integrated LF changes across 
the LF bump as predicted by the models. This 
implies that the H-profile (and thus fuel con- 
sumed in the earlier stages) is given reasonably 
accurately by standard models. 

• The LF is well fit in the SGB region. There is 
no indication of a problem with the canonical 
models as had been earlier found for low metal- 
licity clusters. 

• Below the TO we find evidence for a "flat" IMF 
[x « 0) and evidence for mass segregation which 
is well fit by King models. 



The two best observed GGC LFs, that of M3 pre- 
sented here and M5 by SBSH, seem to be a vindi- 
cation of the quality of theoretical models. In con- 
trast, the anomalous LF of M30 first found by Bolte 
(1994) has inconveniently persisted in later studies 
(Bergbush 1996, Guhathakurta et al. 1998, Sandquist 
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et al. 1999). Why should the theory work so well in 
some cases and fail in others? One can explore possi- 
ble deficiencies in the models as in VandenBerg et al. 
(1998). However, we feel that the ultimate resolution 
of the problem lies in the aquisition of more high qual- 
ity LFs to determine how common anomalous LFs are 
and what cluster parameters they correlate with. 
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6. FIGURE CAPTIONS 



360". 



Fig. 1.— The V, V - I CMD for the HST global 
sample. The sample contains more than 37,000 stars. 



Fig. 13. — Comparison between the X value (com- 
puted as in Bolte 1989) and the dynamical models 
from Pryor et al. 1986. 



Fig. 2.— Map of the WFPC2 field of view The cir- 
cle (at r - 50") delineates the INNER and OUTER 
regions. 

Fig. 3. — The smoothed interpolating curves of com- 
pleteness as a function of V magnitude. 

Fig. 4. — The zoomed LF with and without com- 
pleteness corrections. The Main-Sequence- Turn-Off 
point is labeled. 

Fig. 5. — CMDs of stars selected for the LF in each 
of the three adopted samples: PHOTO, CCD, HST, 
respectively. 

Fig. 6. — Differential Luminosity Function of the 
global sample (PHOTO+CCD+HST). The names we 
use to refer to various parts of the LF are indicated. 

Fig. 7. — Differential luminosity function of the 
global sample compared with theoretical models from 
SCL97 (panel a) and V97 (panel b). The chosen val- 
ues for distance modulus, age and chemical composi- 
tion are indicated. 



Fig. 8. — Integrated luminosity function of the global 
sample compared with theoretical models by SCL97. 



Fig. 9. — The differential LF of the global sample 
is compared with SCL97 theoretical models varying 
metallicity, distance modulus and age. 

Fig. 10. — The LF-bump region of the integrated 
(panel a) and differential (panel b) LF: The position 
of the LF-bump is indicated. The solid line is the ob- 
served LF. The dashed line is based on SCL97 models. 
The short dashed lines in (panel a) are straight line 
fits to the LF above and below the LF-bump and are 
drawn to aid the eye in seeing the change in slope. 

Fig. 11.— The dependence of the SGB LF on 
WIMPs, He diffusion, [a/Fe]. 

Fig. 12. — A comparison of the LFs in three separate 
radial regions: r < 20"; 50" < r < 100"; 210" < r < 
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Mean ridge line for the HST sample. 
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Table 5 

Luminosity function of the global sample (PHOTO+CCD+HST). 
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Table 6 

Chi-square results between observations and SCL models 
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This figure "f5.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9905193vl 
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